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Abstract 

^, 1 In this paper, coordinated beamforming based on relaxed zero forcing (RZF) for K transmitter- 

^ . receiver pair multiple-input single-output (MISO) and multiple-input multiple-output (MIMO) interfer- 

K*" ' ence channels is considered. In the RZF coordinated beamforming, conventional zero-forcing interference 

00 ■ 
lO ' leakage constraints are relaxed so that some predetermined interference leakage to undesired receivers 

i> : 

is allowed in order to increase the beam design space for larger rates than those of the zero-forcing 

CO . (ZF) scheme or to make beam design feasible when ZF is impossible. In the MISO case, it is shown 

O ■ 

PsJ ' that the rate-maximizing beam vector under the RZF framework for a given set of interference leakage 

levels can be obtained by sequential orthogonal projection combining (SOPC). Based on this, exact 

, , , and approximate closed-form solutions are provided in two-user and three-user cases, respectively, and 

;h ■ an efficient beam design algorithm for RZF coordinated beamforming is provided in general cases. 

cd ■ 

Furthermore, the rate control problem under the RZF framework is considered. A centraUzed approach 
and a distributed heuristic approach are proposed to control the position of the designed rate-tuple in 
the achievable rate region. Finally, the RZF framework is extended to MIMO interference channels by 
deriving a new lower bound on the rate of each user 
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I. Introduction 

In current and future cellular networks, handling interference in the network is one of the most critical 
problems. Among the many ways of handling interference, MIMO antenna techniques and base station 
cooperation are considered as the key technologies to the interference problem. Indeed, the 3GPP Long- 
Term Evolution-Advanced considers the base station cooperation and MIMO techniques to mitigate inter- 
cell interference under the name of Coordinated Multipoint (CoMP) lH, |[3l. Mathematically, when each 
mobile station has a single receive antenna and data is not shared among base stations, the system is 
modelled as a MISO interference channel (IC), and extensive research has been conducted on beam design 
for this MISO IC, especially under the assumption of practical linear beamforming treating interference as 
noise. First, Jorswieck et al. investigated the structure of optimal beam vectors achieving Pareto boundary 
points of the achievable rate region of the MISO IC with linear beamforming 0] and showed that any 
Pareto-optimal beam vector at each transmitter is a normalized convex combination of the ZF beam 
vector and matched-filtering (MF) (i.e., maximal ratio transmission) beam vector in the case of two users 
and a linear combination of the channel vectors from the transmitter to all receivers in the general case 
of an arbitrary number of users. The result is extended in []5] to general MISO interference networks 
with arbitrary utihty functions having monotonic property. Moreover, the parameterization for the Pareto- 
optimal beam vector is compressed from K{K — 1) complex numbers Q to K{K — 1) real numbers. In 
addition to these results, other interesting works for MISO ICs include the consideration of imperfect CSI 
|[6ll . shared data Q, second-order cone programming |]8], etc. Although these works provide significant 
theoretical insights into the optimal beam structure and parameterization of Pareto-optimal beam vectors, 
it is not easy to use these results to design an optimal beam vector in the real-world systems, and the 
beam design problem in the general case still remains as a non-trivial problem practically. 

With a sufficient number of transmit antennas, the simplest beam design method for base station 
coordination is ZF, which perfectly eliminates interference leakage to undesired receivers. However, it 
is well known that the ZF method is not optimal in the sense of sum data rate or Pareto-boundary 
achievability, and there have been several ideas to enhance the ZF beam design method. In the case of 
multi-user MISO/MIMO broadcast channels, the regularized channel inversion (RCI) lH and the signal-to- 
leakage -plus-noise (SLNR) method fTO'l were proposed for this purpose. In particular, the SLNR method 
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maximizes the ratio of signal power (to the desired receiver) to leakage (to undesired receivers) plus noise 
power, and its solution is given by solving a generalized eigenvalue problem. The SLNR method can 
easily be adapted to the MISO/MIMO IC. Recently, Zakhour and Gesbert rediscovered this method in 
the context of MISO IC under the name of the virtual signal-to-interference-plus-noise (SINR) method, 
and have further (and more importantly) shown that this method can achieve any point on the Pareto 
boundary theoretically, but practically can achieve one uncontrolled point on the Pareto boundary of the 
achievable rate region in the case of twofj users |[T2l . ||T31 . 

Another way of generalizing ZF in MISO IC was proposed by relaxing the ZF leakage constraints to 
undesired users in ||T41 . lITSl . [[T]. First, Shang et al. showed that all boundary points of the achievable rate 
region of MISO IC with single-user decoding can be obtained by linear beamforming ||T41 . by converting 
the non-convex weighted sum rate maximizing precoder design problem into a set of separate convex 
problems by taking a lower bound on the achievable rate of each user under the relaxed ZF (RZF) 
framework. This method was further investigated by Zhang and Cui |[T5l . who showed that separate 
rate optimization under the RZF framework with a set of well-chosen interference leakage levels to 
undesired users is Pareto-optimal for MISO ICs in addition to being sum-rate optimal. In |[T], Lee et 
al. extended the RZF framework to the case of MIMO IC. In this RZF beamforming framework, each 
transmitter maximizes its own rate under interference leakage constraints to undesired receivers. The idea 
is based on the simple observation that the ZF beam design method overreacts to inter-cell interference 
by completely nulling out the interference. Most receivers (i.e., mobile stations) that are affected by inter- 
cell interference are cell-edge users, and thus, thermal noise remains even if the inter-cell interference is 
completely removed. Thus, it is unnecessary to completely eliminate the inter-cell interference and it is 
sufficient to limit the inter-cell interference to a certain level comparable to that of the thermal noise. By 
relaxing ZF interference constraints, we do not need the condition that the number of transmit antenna 
is larger than or equal to that of receivers and have a larger feasible set yielding a larger rate than that of 
the ZF scheme. In this paper, we explore and develop this RZF idea fully in several aspects to provide a 
useful design paradigm for coordinated beamforming (CB) for current and future cellular networks. The 
contributions of the paper is summarized as follows: 

• In the MISO IC case, a new structural representation of optimal beam vector for RZF coordinate 
beamforming is derived. 

'it can be shown that the virtual SINR (or SLNR) method can theoretically achieve any Pareto-optimal point in the general 
MISO IC case, too. See the appendix of (TTI . 
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• In the MISO IC case, based on the new structural representation, the sequential orthogonal projection 
combining (SOPC) method for the RZF beam design is proposed. In the case of K = 3, an approximate 
closed-form solution is provided. 

• In the RZF framework, the allowed interference leakage levels to undesired receivers at each transmitter 
are design parameters, and the rate-tuple is controlled by controlling these interference leakage levels. A 
centrahzed algorithm and a fully distributed heuristic algorithm are provided to control the location of the 
designed rate-tuple (roughly) along the Pareto boundary of the achievable rate region. The controllability 
of rate is a desirable feature in network operation since the required data rate of each transmitter-receiver 
pair may be different from those of others in practice, as in an example that one user is a voice user and 
the others are high rate data users. 

• Finally, the RZF CB (RZFCB) is extended to the MIMO IC case. In the MIMO case, a new lower 
bound on each user's rate is derived to decompose the beam design problem into separate problems at 
different transmitters, and the projected gradient method lfT6ll is adopted to solve the MIMO RZFCB 
problem. 

Notations and Organization In this paper, we will make use of standard notational conventions. Vectors 
and matrices are written in boldface with matrices in capitals. All vectors are column vectors. For a matrix 
A, A^, ||A||, ||A||i?, tr(A), and |A| indicate the Hermitian transpose, 2-norm, Frobenius norm, trace, 
and determinant of A, respectively, and C(A) denotes the column space of A. I„ stands for the identity 
matrix of size n (the subscript is omitted when unnecessary). IIa = A(A^A)~^A^ represents the 
orthogonal projection onto C(A) and 11^ = I — IIa. For matrices A and B, A > B means that A — B 
is positive semi-definite, [ai,--- ,3.1] or [a^]^^ denotes the matrix composed of vectors ai,--- ,ai. 
X ~ CM{ti, XI) means that x is circular-symmetric complex Gaussian-distributed with mean vector /j, 
and covariance matrix Xl. M, M+, and C denote the sets of real numbers, non-negative real numbers, and 
complex numbers, respectively. For a set A, \A\ represents the cardinality of the set. 

The remainder of this paper is organized as follows. The system model and the preliminaries are 
provided in Section JIl In Section |llll the RZFCB in MISO ICs is formulated, and its solution structure 
and a fast algorithm for RZFCB are provided. In Section JV] the rate-tuple control problem under the 
RZFCB framework is considered and two approaches are proposed to control the designed rate-tuple. 
The RZFCB problem in MIMO ICs is considered in Section |Vl followed by conclusions in Section IVll 
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II. System Model and Preliminaries 

In this paper, we consider a multi-user interference channel with K transmitter-receiver pairs. In the 
first part of the paper, we restrict ourselves to the case that the transmitters are equipped with N antennas 
and each receiver is equipped with one receive antenna only. In this case, the received signal at receiver 
i is given by 

K 
Vi = ^fi^iSi + ^ ^ij^jSj + rii, (1) 

where hij denotes the A^ x 1 (conjugated) channel vector from transmitter j to receiver i, and v^ and 
Sj are the A^ x 1 beamforming vector and the scalar transmit symbol at transmitter j, respectively. We 
assume that the transmit symbols are from a Gaussian code book with unit variance, the additive noise 
Ui is from CM{0, of), and each transmitter has a transmit power constraint, ||vj|p < Pj, i = 1, • • • ,K. 
The first term on the right-hand side (RHS) of ([T]) is the desired signal and the second term represents 
the sum of interference from K — \ undesired transmitters. Under single-user decoding at each receiver 
treating interference as noise, for a given set of beamforming vectors {vi,--- ,vj^} and a channel 
realization {hjj}, the rate of receiver i is given by 

i?.(vi, • • • , V;,) = log (l + ^.,^"";[h^.p ) • (2) 

Then, for the given channel realization, the achievable rate region of the MISO IC with transmit beam- 
forming and single-user decoding is defined as the union of the rate-tuples that can be achieved by all 
possible combinations of beamforming vectors under the power constraints: 

7^:= U (i?i(vi,--- ,vk), ••• , i?i^(vi,--- ,vi^)). (3) 

/ v,:v.6C", 1 

ll|v.||^<P,, ^<i<K] 

The outer boundary of the rate region Tt is called the Pareto boundary of IZ and it consists of the rate- 
tuples for which the rate of any one user cannot be increased without decreasing the rate of at least one 
other user |@]. 

At each transmitter, the interference to undesired receivers can be eliminated completely by ZF CB 
(ZFCB). Due to its simplicity and fully distributed nature, there has been extensive research on ZFCB, 
e.g., llT7l - |[T9l . The best ZF beamforming vector at transmitter i can be obtained by solving the following 
optimization problem: 

V* = arg max log 1 + ' " '' (4) 

subject to |h^Vj| =0, V j / i and ||vi|p < Pi. 
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Here, |h^Vj| = is the ZF leakage constraint at transmitter i for receiver j. If A^ > K, the problem 
dH) has a non-trivial solution and the solution is given by vf ^ = clli . . u i hjj for some 

scalar c satisfying the transmit power constraint. In this paper, however, we do not assume that N > K 
necessarily as in the ZF beamforming, but assume that 

(AA) In the case of N > K, {h.ji,j = !,••• ,K} are linearly independent for each i. In the case 
of A^ < K, the element vectors of any subset of {hj,, j = 1, • • • ,K} with cardinality N are linearly 
independent for each i. 
Assumption (A.1) is almost surely satisfied for randomly realized channel vectors. 

III. RZF Coordinated Beamforming in MISO Interference Channels 

A. Formulation 

Although the ZFCB provides an effective way to handling inter-cell interference, the ZFCB is not 
optimal from the perspective of Pareto optimality, i.e., the rate tuples achieved by ZFCB are in the 
interior of the achievable rate region 11201 . and requires the condition A^ > K. As mentioned before, 
even with such complete interference nulling, there exists thermal noise at each receiver, and thus, a 
certain level of interference leakage comparable to the power of thermal noise can be allowed for better 
performance. In the MISO IC case, the RZF leakage constraint at transmitter i for receiver j is formulated 
as follows: 

Ihjfvjp < ajicrj, Vz, j / i, (5) 

where Uji > is a constant that controls the allowed level of interference leakage from transmitter 
i to receiver j relative to the thermal noise level (t^ at receiver j. When aji = for all j ^ i, the 
RZF constraints reduce to the conventional ZF constraints. When aji > 0, on the other hand, the ZF 
constraints are relaxed to yield a larger feasible set for Vj than that associated with the ZF constraints 
and due to this relaxation the condition A^ > i^ is not necessary anymore. 

Under the RZF framework, the power of interference from unde sired transmitters at receiver i is upper 
bounded as 

Ef=l,j^^ |hf- vjf < Zj^i ^ijaf =: eiul (6) 

^In the RZF scheme, {aji ,j,i = l,---,K,j^i} are system design parameters that should be designed properly for optimal 
performance. The practical significance of the parameterization in terms of the interference leakage levels will be clear in Section 
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Therefore, a lower bound on the rate of user i under RZF is obtained by using ^ as 

The lower bound on the rate at each receiver does not depend on the beamforming vectors of undesired 
transmitters and thus, exploiting the RZF constraints, we can convert the intertwined coordinated beam 
design problem into a set of separate problems for different users based on the lower bound ril4i . The 
separate problem for each transmitter based on RZF is given as follows |[T4l . |[T5l : 
Problem 1: For each transmitter i G {1, • • • , K}, 

maximize log 1 + -j^-^ — ^^-^ ) (8) 



(1 + e.)< 
subject to |h^Vip < ajio], Vj / i, (9) 

\\^^f<Pi■ (10) 

Then, due to the monotonicity of the logarithm. Problem [U is equivalent to the following problem: 
Problem 2 {The MISO RZFCB problem): For each transmitter i G {1, • • • , K}, 

maximize |hjjVjp (11) 

Vi 

subject to |h^Vjp < ajia'j, Vj / i, (12) 

||v.,f<Pi. (13) 

From now on, we will consider Problem |2] (the RZFCB problem) and refer to the solution to Problem |2] 
as the RZF beamforming vector. 

B. The Optimality and Solution Structure of RZFCB in MISO Interference Channels 

In this subsection, we will investigate the optimality and structure of the solution to Problem |2] We 
start with the optimality of the RZFCB scheme. Without inter-cell interference, it is optimal for the 
transmitter to use the MF beam vector with full transmit power. However, with inter-cell interference, 
such a selfish strategy leads to poor performance due to large mutual interference pOl. Thus, to enhance 
the overall rate performance in the network, the beamforming vector should be designed to be as close 
as possible to the MF beam vector without giving too much interference to undesired receivers, and this 
strategy is the RZFCB in Problem |2] (or Problem [T]equivalently). The optimality of the RZFCB is given 
in the following theorem of Shang et al. |[T4l or Zhang and Cui |[T5l . 
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Theorem 1: |[T5l Any rate-tuple {Ri, • • • , Rk) on the Pareto boundary of the achievable rate region 
defined in Q can be achieved by the RZFCB if the levels {aijaf, Vi, j / i} of interference leakage are 
properly chosen. 

Proof: See Proposition 3.2 in fT5l . ■ 

Surprisingly, the separate beam design based on the rate lower bound in Problem |2] can achieve any 
Pareto-optimal point of the achievable rate region if the interference relaxation parameters are well 
choseno It was also shown that Problem 2 and the approach in |'5l| are two different approaches to the 
same multi-objective optimization problem ETI . Due to Theorem [T] in the MISO IC case, the remaining 
problems for the RZFCB are i) to construct an efficient algorithm to solve the RZFCB problem for given 
{aijaf, yi,j 7^ i} and ii) to devise a method to design {aijaf, 'ii,j 7^ i} for controlling the location 
of the rate-tuple along the Pareto boundary of the achievable rate region. We will consider Problem |2] 
for given {aijaf, \/i,j / i} here and will consider the rate control problem in the next section. 

First, we will derive an efficient algorithm for obtaining a good approximate solution to Problem |2] for 
given {aijaf, \/i,j^i}.Todo this, we need to investigate the solution structure of the RZFCB problem. 
Instead of solving Problem [T] as in lITSl (this becomes complicated due to the logarithm), we here solve 
Problem |2j which is equivalent to Problem [T] Note that Problem |2] is not a convex optimization problem 
since it maximizes a convex cost function under convex constraint sets instead of minimizing the cost. 
However, Problem |2] can be made an equivalent convex problem by exploiting the phase ambiguity of 
the solution to Problem |2] and making h^Vj real and nonnegative without affecting the value of |h^Vj| 
as follows II22I : 

Problem 3: For each transmitter i G {1, • • • , K}, 

maximize hjjVj (14) 

subject to |hjjvjp < ajia^, Vj / i, (15) 

h^f<P^, (16) 

hf V, > 0. (17) 

Here, the constraint dTTT i implies imag(h^Vj) = and due to this constraint, maximizing |h^Vjp is 
equivalent to maximizing h^Vj. 

^The beamforming vectors from Problem |2] are necessary to achieve any point on the Pareto boundary but not sufficient. Not 
any choice of parameters {oij } leads to a point on the Pareto boundary. 
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Lemma 1: Let v°^ be a solution of the RZFCB problem (i.e., Problem |2]| for transmitter i. Then, v^^ 
is represented as follows: 

v°^* = Ciih.ii + ^ Cjihji (18) 

for some {cji G C : j G T^ U {i}}, where Ti := {j : |h|[vf *|2 = Ujia^}, ||vf *f = Pi for iV > i^, and 
||vf *f < -Ri for iV < /s:. 

Proof: Proof is based on the equivalent formulation in Problem |3] Since Problem |3] is a convex 
optimization problem, the optimal solution can be obtained by the Karush-Kuhn-Tucker (KKT) conditions. 
The Lagrangian of Problem |3] for transmitter i is given by 

£(vi,A,/i,i/) = -hfv, (19) 

K 

+ Y. A,(|hjjv,|2 - a,ia]) + /i(||v,f - P,) - i^hf v,, 

where A := {Aj > : j = 1, • • • , i — 1, i + 1, • • • , K} and fi,u > are real dual variables. With optimal 
dual variables A*, ^*, and u*, the (complex) gradient of the Lagrangian should be zero at v°^ , i.e., 

= Vv*/:(v„ A^;u^l/'^)|^^^^„,. (20 

where Fj := {j : A^ > 0} and Vv* is the conjugate Wirtinger gradient. From the complementary 
slackness condition, A^ > only when |h^Vjp = (Xjiaf. Also, from the complementary slackness, we 
have u* = 0. Otherwise, h^v°^ = and thus no rate is provided to user i. Thus, the gradient of the 
Lagrangian becomes zero if and only if 

/l + E,er.A*h,,h|f)vf. (21) 

If Q := (/_i*I + X^ver Athjjh^) is singular, then v°^ exists if and only if ha G C(Q). However, the 
condition ha G C(Q) does not occur almost surely for randomly realized channel vectors, which is 
assumed here. Therefore, Q should have full rank for the existence of v°^ and the corresponding v°^ 
has two different forms according to the optimal dual variable fi*. 

i) fi* > 0: This corresponds to the case in which the transmitter uses full power, i.e., \\v°^ |p = Pi. 
In this case, the optimal solution is given by 

vf = (/.-I + Zjer. \^h,,hjj) "' h,,. (22) 
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By applying the matrix inversion lemma recursively, it can be shown that v°^ is a linear combination 
of {hji : J G r^ := Fj U {i}}. Thus, the solution is represented as (ITSl l. 

a) fi* = 0: This case corresponds to the case in which full power is not used at transmitter i. In this 
case, Q = X^igr ^l^ji^fi- "^^^ matrix Q in this case is non-singular if and only if |rj| > A^ (i.e., 
K > N) under the assumption (A.l), and the corresponding solution is given by 

vr=(E,er.A,*h,,h^^)"'h,,. (23) 

In this case, {hjj, j G Fj} alone span C^ fully and it is therefore clear that the solution is represented 
as (ITSl l. Indeed, any subset of {hji,j = !,••• ,K} with cardinality N forms a full basis for C^ under 
the assumption (A.l) in this case. 

Furthermore, when N > K, vf ^ is feasible and thus, we can always increase power and rate without 
causing interference to the undesired receivers. Therefore, the optimal solution uses full power, i.e., 
||v°^*|p = Pi when N > K. On the other hand, when N < K, we can have either n* > (||v°^*|p = Pi) 
or^* = 0(||vf*||2<P,). ■ 

The solution to RZFCB for a given set of interference relaxation levels is a linear combination of the 
desired channel and a subset of interference channels for which the RZF constraint (IT2l ) is satisfied with 
equality. Furthermore, it was shown that the interference leakage levels should be designed to make the 
RZF interference leakage constraints be satisfied tightly in order to achieve a point on the Pareto boundary 
lITSl . In this case, Fj = {1, • • • , K}\{i} and thus, the RZF beam structure in Lemma d] coincides with the 
Pareto-optimal beam structure derived by Jorswieck et al. in lH. Now, based on Lemma [T] we present a 
new useful representation of v°^ that provides a clear insight into the RZFCB solution and a basis for 
fast algorithm construction. 

Theorem 2: For transmitter i, the RZFCB solution can also be expressed as 

opt ^^^ ^ njh^^ "A„^hn 

"«« ■'-'■A,'-^** ■'-'■A- '■'■iiW 

where Cj G C, j = 0, 1, • • • , |rj| and Aj is constructed recursively as 

Aj := [Aj-_i, hp_(^.)^.], j = 1, • • • , |fi|. (25) 

Here for convenience we let Aq be an A^ x 'matrix'. Fj is a set made by permuting the elements of 
Fj according to an arbitrary order, and Fj(j) denotes the j'-th element of Fj. 

Proof: From Lemma[T] we know that v°^ G C([hjj]jgr')- Proof of the theorem is given by showing 
the equivalence of the two subspaces C{\h.ji]j^Y[) and C([hjj,n;^ hjj, • • • , H;^^ hji]). 
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Case (i). |rj|(= |rj|) < A^ — 1.' In this case, {hjj, 11;^ hjj, • • • , 11;^ _ hjj} are linearly independent. 
This is easily shown by replacing 11^ with I — IIaj and by using the linear independence of {hjjjjgr'- 
Thus, the dimension of C([hjj,n;^ ha,-- , 11;^^ hjj]) is (|rj| + 1), which is the same as that of 
C([hjj]jgrj)- Now, consider the projection of any vector in C([hjj,n;^ hjj,--- ,11;^^ hjj]) onto the 
orthogonal complement of C([hjj]jgr')- 

n[k,,],,r, (cohn + cini^hjj + - - - + C|p^|ni|^^ hjj) 

= nfi^^.^]^^^, (cohii + ci(I-nAi)hjjH hC|p^|(I- nA|f;^|)hii) 

= ni I (cohji + ci(I - Hhp )hjj + - - - + C|p,|(I - n ,|F,i)hii) 



= 0. (26) 

By (|26l) the orthogonal complement of C( [h;j. Uj^ hjj,--- ,11;^^ hjj]) is included in that of C([hjj]jgpj), 
but C( [hjj, 11;^ hjj, - - - , 11;^ ^ hjj] ) and C([hjj]jgr;) have the same dimensions. Thus, the two orthogo- 
nal complements are the same, and hence, the two subspaces themselves are the same. Consequently, for 
any Cjjhjj + J^jer '^ji^ji ^^^^ arbitrary {cji G C : j G F'-}, there exists some {cj G C : < j < jFjl} 
s.t. 

opt u , X^ u hjj ^^ ^^A,"** 

V/ = Cjjhjj + 2^ Cjihji = CQ——r + 2^ Cj . 

Case (ii). \Ti\ > N: In this case, both {h^j, j G Fj} and {hjj, 11^ hjj, - - - , 11^ ^ hjj} span the whole 
C^. Thus, the claim is trivially satisfied. ■ 

Theorem |2] states that the RZF solution is a linear combination of vectors that are obtained by projecting 
the desired channel vector onto the orthogonal complements of a series of subspaces spanned by the 
channels from the transmitter to the undesired receivers. Furthermore, the series of subspaces are obtained 
by sequentially including one additional interference channel vector at a time, as shown in (|25T l. Soon, it 
will be shown that, to obtain the RZF solution to Problem [2l the order of interference channel inclusion 
for constructing AjS in Theorem |2] is determined by the set of allowed interference levels and the channel 
realization. 

C The Sequential Orthogonal Projection Combining Method and Closed-Form Solutions 

In this subsection, we propose an efficient beam design method for RZFCB that successively allocates 
the transmit power to certain vectors obtained by sequential orthogonal projection of the desired channel 
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vector onto monotonically decreasing subspaces. Furthermore, we provide the closed-form solution to the 
RZFCB problem in the two-user case and an approximate closed-form solution in the three-user case. 

To obtain the RZF beamforming vector under given interference relaxation constraints for a given 
channel realization, Problem |2] should be solved. One can use a numerical method ll23l . as in lITSl . 
However, such a method requires a numerical search for determining the Lagrange dual variables sat- 
isfying the RZF constraints and the transmit power constraint. To circumvent such difficulty and to 
increase the practicality of the RZFCB, we exploit Theorem |2] to construct an efficient method to find the 
RZFCB solution. Theorem |2] provides us with a very convenient way of obtaining the RZFCB solution 
for given interference leakage levels for a given channel realization; we only need to find Fj and complex 
coefficients {q} in (l24l ) for each transmitter. The idea is based on the fact that the RZF beamforming 
vector should be designed to be as close as possible to the MF beam vector under the interference 
leakage constraints for the maximum rate under RZF, as described in Problem |2l Hereafter, we will 
explain how the coefficients {q} and the matrices {Aj} in Theorem |2] can be obtained to maximize the 
rate under the RZF interference and power constraints. Consider transmitter i without loss of generality. 
For the given transmit power constraint ||vj||2 < Pi, it may not be possible to allocate all of the transmit 
power to the MF direction ha because this allocation may violate the RZF constraints. The rate greedy 
approach under the RZF constraints for a given channel realization is explained as follows. First, we 
should start to allocate the transmit power to the direction of ha by increasing cq with some phase 
until this allocation hits one of the RZF constraints with equality, i.e., the interference level to one of 
the undesired receivers reaches the allowed maximum exactly. (In the case that the allowed interference 
levels to all undesired receivers are the same, this receiver is the receiver whose channel vector has the 
maximum inner product with ha.) The index of this receiver is rj(l). At this point, transmitter i cannot 
allocate the transmit power to the direction ha anymore since this would violate the RZF constraint for 
receiver rj(l). Since the RZF constraints for other undesired receivers are still met with strict inequality, 
transmitter i can still cause interference to the remaining receivers. Thus, for the maximum rate under 
the RZF constraints, transmitter i should now start to allocate the remaining power to the direction of 
n^^hjj, where Ai = [hp ^^ J, until this allocation hits another RZF constraint with equality. The index 
of this receiver is rj(2). (Note that 11^ ha is the direction of maximizing the data rate without causing 
additional interference to receiver rj(l).) Now, transmitter i cannot cause interference to receiver ri(2) 
in addition to receiver rj(l) anymore. Therefore, at this point, transmitter i should start to allocate its 
remaining power to the next greedy direction 11^ ha, where A2 = [hf /^^n ■, hp .2-1 J- This greedy power 
allocation without violating the RZF constraints should be done until either all the transmit power is 
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used up (fi* > in Lemma [D or we cannot find a new direction that does not cause interference to 
the users that are already in the set Fj (/i* = in Lemma [T]). When N > K and transmit power still 
remains even after hitting all the K — 1 interference leakage constraints with equality, from then on, 
all the remaining power should be allocated to the ZF direction. Fhis coincides with our intuition that 
ZF is optimal at a high signal-to-noise ratio (SNR) in the case of A^ > K. On the other hand, when 
all the transmit power is used up before reaching the remaining interference constraints with equality, 
the corresponding remaining interference channel vectors do not appear in the solution. The final RZF 
solution is the sum of these component vectors and has the form in (l24l) . In this way, the RZFCB 
solution can be obtained by combining the sequential projections of the desired channel vector ha onto 
the orthogonal complements of the subspaces C(Ai) C ••• C C(A|p ,). Thus, we refer to this beam 
design method as the sequential orthogonal projection combining (SOPC) methodu By Theorem [T] the 
SOPC strategy with a well chosen set of interference relaxation levels is Pareto-optimal for MISO K-pair 
interference channels with single-user decoding. 

An interesting interpretation of the SOPC strategy is in an analogy with the water-filling strategy. 
The water-filling strategy distributes power to resource bins according to the effectiveness of each bin, 
and the power fills into the bin with the lowest noise level (or the most effective bin) first. Similarly, 
the SOPC strategy allocates power to the most effective direction first and then the next most effective 
direction when the first direction cannot accommodate power anymore. This procedure continues until 
either the procedure uses up the power or it cannot find a new feasible direction. So, the SOPC strategy 
can be viewed graphically as pouring water on top of a multi-tiered fountain, as illustrated in Fig. [T] 
The relationship of the RZFCB/SOPC design and the two-user result by Jorswieck et al. llH is explained 
in Fig. |2] In the two user case, Jorswieck et al. have shown that a Pareto-optimal beam vector is a 
convex combination of the MF beam v^^^ and the ZF beam vf ^ satisfying the power constraint, i.e.. 



Vj = y/Pi ||/yMi. ' TQl;^'\yijj. ' || , where < Aj < 1. Thus, the feasible set of optimal beam vectors is the 
arc denoted by F in Fig. |2] All the points on this arc can be represented by the sum of the two vectors 
in red, and the size of the component vector in the MF direction is determined by its projection onto 

''The rate optimality of the SOPC strategy under the RZF constraints is straightforward to see. Suppose that we are given any 
beam vector that is a linear combination of {hji}, satisfies the RZF interference and power constraints but is not the SOPC 
solution. Then, the vector can still be represented in terms of the SOPC basis in Theorem |2] and some of the basis component 
vectors with larger iimer product with the MF direction do not satisfy the RZF constraints with equality. Thus, the rate can be 
increased by allocating power from the basis component vector with smaller inner product with the MF direction to the basis 
component vector with larger inner product with the MF direction until the RZF constraints are satisfied with equality. 
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Fig. 1. The SOPC strategy in tlie case of N > K: Water-pouring on a multi-tiered fountain. 
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Fig. 2. SOPC interpretation of the two-user result. 

C{h.2i), i.e., the allowed interference level to the other receiver in the RZF context. Thus, the two-user 
result by Jorswieck et al. can be viewed as a special case of the SOPC strategy when the number of 
users is two. The key difference is the parameterization; ai2 and 021 are the parameters in the RZF 
framework whereas the linear combining coefficients Ai and A2 are the parameters in |@]. 

Now consider the detailed implementation of the SOPC method. Before considering the general case 
of an arbitrary number K of users, we consider simple two-user and three-user cases. Here, we restrict 
the combining coefficients {q} to the set of real numbers. It will shortly be shown that the performance 
loss caused by restricting {q} to real numbers is negligible. Furthermore, it is the optimal solution of 
the RZFCB when K = 2. For simplicity, we only provide the solution for transmitter 1. The solutions 
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for other transmitters can be obtained in a similar way. 

Proposition 1: The closed-form SOPC solution in the two-pair MISO IC case is given by 



Vl 



/JYvf ^, if P < ^^^i^ 



"21^1 

CoVi^^ + '^ivf^, otherwise, 



where vf ^ = r^, vf ^ = ^^^7^, ^o = \fJ^, and 6 = -p^o + ^W^W^^^T)- Here, 

I!-'--'-!! II hoill V'211 I 

(vf^)^vf^ = ||ni^ hii||/||hii||GM+. 



P 



Proof: Proof of Proposition [T] can be found in lITTI . ■ 

Now, we consider the case of K = 3. This case is particularly important when the hexagonal cell structure 
is used and three cells are coordinating their beam vectors. In the case of K = 3, the solution can have six 
different forms depending on the transmit power and channel realization. We will provide the closed-form 
solution under the real coefficient restriction for transmitter 1 in the case that the interference leakage 
to receiver 3 reaches the allowed level before the interference leakage to receiver 2 reaches the allowed 
level. (For this, we should first take inner products (h2i,hii) and (h3i,hii) and compare the ratio of 
their magnitudes with some threshold. The solutions of the other case and of other users can be derived 
in the same manner.) 

Proposition 2: For K = 3 and . -j/ ij^L > ^^ i , the closed-form SOPC solution with the restriction to 

irijiVi p a2iO"2 

/ 2 

real coefficients at transmitter 1 is given in (l27l ). In (l27l) . /3o^ /3i, /3i and /32 are given by /3o = ^tttSt. 

/3i = -af3o + ^Pi - (1 - a^)i3l I3[ = i(-d/3o + ^/d^M^^^ibi3i^^^^) , and /32 = -(//3o + 



e(3[) + V(//?o + e/3i)2 - (2a/3o/3; + /Sq" + /3f - Pi), where a := Re{(vf ^, ni^^^hii/||ni^^^hii||)} 
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Proof: Proof of Proposition |2] can be found in lITTI . ■ 

In the case of K > 3, it is cumbersome to distinguish all possible scenarios for deriving an explicit 
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Number of users 

Fig. 3. Average sum rates of the exact RZFCB solution and the proposed SOPC algorithm with real coefficients. (Here, N = K 
and the average sum rate is obtained over 50 i.i.d. channel realizations.) 



SOPC solution. Thus, we propose an algorithm implementing the SOPC strategy with real combining 
coefficients in Table Jl In the general case of i^ > 3, the implementation of the SOPC algorithm can be 
simplified by the known result in the Kalman filtering theory, provided in the following lemma. 

Lemma 2 (Sequential orthogonal projection ^24^): Let ^ be a Hilbert space with norm || • || and inner 
product (•,•). Consider x G H and a closed linear subspace Aj of T-L. For some y G H but y ^ Aj, the 
following equality holds 



Ha. + iX = n[A„y]X 



(28) 



rr (x-nA,x, y-nA,y) , 

nA,x + n — T^ — f?5 (y - nA,y). 



l|y-nA,yP 

Since we need to compute 11;^ hjj = (I — nAj)hjj in the SOPC algorithm. Lemma |2] can be applied 
recursively by exploiting the fact Aj = [Aj_i,hp. .^ J. Thus, we only need to compute nAj_ihp, n^ 
for each j G {1, 2, • • • ,K}. The proposed algorithm in Table U computes the direction and size of the 
component vector for SOPC directly in each step. 

The proposed SOPC solution based on real coefficients is a sub-optimal solution to the RZFCB problem 
in the case of K > 3. However, the performance loss between the optimal RZFCB (or exact SOPC) 
beamforming vector and the proposed SOPC solution based on real coefficients is insignificant for a 
wide range of meaningful SNR values, as seen in Fig. |3] Thus, practically, the proposed SOPC solution 
can be used with negligible performance loss. Note that the necessary computations for the proposed 
SOPC solution are a few inner product and square root operations and the complexity of the SOPC 
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TABLE I 

The sequential orthogonal projection combining algorithm. 



Given channel realization 
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,i,j = I,--- 


,K}, 
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{aj,a] : i,j = I,--- ,K, 
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i}, and maximum transmit 
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I,... 
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2. ^p is a positive solution of 


■ ||v, + /ipu||i = 


= P^, 


i.e., 









Aip := -pp + \JpI - (||vi||i - P,) 



where pp — _Re(u^Vi). 

3. /ij is a positive solution of |hj^(vi + /iju)p = otjicr'^ 
for each j G $i, i.e., 

''j -^ |hHu|2 

where Pj — -Re(vf^hjihjfu). 

4. Obtain /i^ = min{/ij} and j* = argmin{/^j}. 

5. If ^p > p*, Vi = Vi + p*vi, A = [A.hj.i], $i = $A{J*}' 
k = k + 1, and go to step 1. 

If Pp < Pj , "Vi = Wi -\- ppU. Terminate iteration. 
end 



method is simply 0{N), where A'^ is the number of transmit antennas at the transmitter. The proposed 
SOPC method reduces computational complexity to obtain an RZF solution by order of hundreds when 
compared to the ellipsoid method for the RZFCB solution used in ifTSl . as shown in Fig. IH and the 
solution procedure can easily be programmed in a real hardware. 

IV. Rate-tuple Control 

In the previous section, we provided an 0(A^)-complexity algorithm to solve the RZFCB problem for a 
given set {aji} of interference relaxation parameters. Now, we consider how to design these parameters. 
We first provide a centralized approach to determine {ojj} with the aim of controlling the rate-tuple 
along the Pareto boundary of the achievable rate region and then a fully-distributed heuristic approach 
that exploits the parameterization in terms of interference relaxation levels in RZFCB and is able to 
control the rate-tuple location roughly along the Pareto boundary of the achievable rate region. 
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Fig. 4. Computational complexity for RZFCB beam design: Ellipsoid method 1151 . 1231 versus SOPC (A'^ — K and SNR=5 
dB) 



A. A Centralized Approach 

By Theorem [U with a set of well chosen allowed interference leakage levels, the RZFCB can achieve 
any Pareto-optimal point of the rate region. However, the problem of designing the interference leakage 
levels {uji] in the network remains. Under the RZFCB framework, in ifTSi, a necessary condition for the 
interference relaxation parameters at each receiver to achieve a Pareto-optimal point was derived. Based 
on the necessary condition, the authors proposed an iterative algorithm that updates the interference 
relaxation parameters. Although the algorithm in lITSi is applicable to general A'-user MISO interference 
channels, it cannot control the rate-tuple location on the Pareto boundary to which the algorithm converges. 
To control the rate-tuple to an arbitrary point along the Pareto-boundary of the achievable rate region, we 
here apply the utility function based approach in ||251 to the RZF parameterization in terms of interference 
leakage levels. Exploiting the fact that the RZFCB can achieve any Pareto-boundary point by adjusting 
{dji}, we convert the problem of finding a desired point on the Pareto boundary of the achievable rate 
region into that of finding an optimal point of the following optimization problem: 



subject to |h^Vip < ajifT|, Vi, j ^ i, 



max 



(29) 



<Pr. 



Vz, 



where u{Ri, • • • , Rk) is the desired utility function and several examples include the weighted sum rate 

u{Ri, • • • , Rk) = Yl WiRi, where Wi> and ^ Wj = 1, the Nash bargaining point u{Ri, • • • , Rk) = 



Hill (Ri-Rf^):. where Rf^ = log2 ( 1 + 2 ,^"\gL/F|2 ) , and the egaUtarian point u{Ri , • • • , Rk) 
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min(i?i, • • • , Rk) l[25l . The optimization ( [29l ) can be solved by an alternating optimization technique. 
That is, we fix all other ajj's except one interference relaxation parameter and update the unfixed 
parameter so that the utility function is maximized. After this update, the next Uji is picked for update. 
This procedure continues until converges. The proposed algorithm is described in detail in Table JI] For 
a given utility function u{Ri, ■ ■ ■ ,Rk), the RZF beam vectors {v^^^} can be obtained as functions of 
{ttji} by the SOPC method, the rate-tuple can be computed as a function of {vf^^} by ©, and finally 
the utility function value can be computed as a function of (i?i, • • • , Rr)- Thus, the utility value as a 
function of {aji} can be computed very efficiently by the SOPC method for the proposed centraUzed 
algorithm, and this fact makes it easy to apply a numerical optimization method such as the interior point 
method to the per-iteration optimization in Table Ull 
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Fig. 5. Convergence of the proposed centralized approach (K = N = 2, Pi = at = 1 foi i = 1, 2). 

Due to the non-convexity of utility functions w.r.t. {oji}, the convergence of the proposed algorithm 
to the global optimum is not guaranteed, but the proposed algorithm converges to a locally optimal 
point by the monotone convergence theorem since the utility function is upper bounded and the proposed 
algorithm yields a monotonically increasing sequence of utility function values. Furthermore, the proposed 
algorithm is also stable by the monotone convergence theorem. Fig. [5] shows the convergence behavior 
of the proposed utility function based algorithm for 10 different channel realizations when K = N = 2, 
Pi = ai = 1 {i = 1,2), and u{Ri,R2) = 2Ri + i?2- It is seen in the figure that the algorithm 
converges in a few iterations in most cases. Fig. |6] shows the convergence behavior of several known 
rate control algorithms for the same setting as in Fig. |5] for one channel realization. The considered 
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TABLE 11 

A CENTRALIZED ALGORITHM FOR DETERMINING {Oji}. 

For given channel realization {hji,i,j = 1, • ■ • ,K}, noise power {cr^,i = 1, ■ • ■ ,K}, 
and a utility function u({i?,}), perform the following procedure to determine interference 
leakage levels {otji}. 

Initialization: {aj; =0, i,j = l,---,K,j^ i}, {R'f =0, i = 1, ■ ■ ■ , K}, e > 0, and 

1 = 1. 

while \u{{a'^,}) - ui{a\T^})\ > e 
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three algorithms converge to the same value eventually in this case. It is also seen in Figs. |7] (a) and (b) 
that the proposed centralized algorithm yields desired points on the Pareto boundary although it is not 
theoretically guaranteed. 



CD 


fl) 


cn 


(/) 




3 


E 


"ri 


3 


f- 


(0 


c- 


t: 


CO 


(11 


J- 


.H- 


n 


J- 

m 




cu 


(0 



1.5 



5^0.5 





- 


1 


Proposed algorithm 

Utschick and Brehmer [1 2-SP] 
Qiuetal[11-SP] 





10" 



10' 
Number of iterations 



10' 
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B. A Distributed Heuristic Approach and Practical Considerations 

The proposed centralized algorithm in the previous subsection requires central processing with the 
knowledge of all {h^j : i,j = 1, 2, • • • , K} and {cr^ : j = 1, 2, • • • , K}. This reduces the practicality 
of the centralized approach when communication among the base stations is limited or experiences 
large delay as in real systems. Note that the RZFCB framework in Problem |2] itself is distributed. 
Transmitter i only needs to know {h.ji,j = 1,2, ••• ,K} and {a'j,j = 1,2, ••• ,K} and needs to 
control {an, • • • , aj-i,j, aiJ^i^i,aKi}- In the RZF framework, heuristic rate control is possible with the 
knowledge of {hji, j = 1, 2, • • • , K} and {a? j = 1, 2, • • • , K} at transmitter i. For fully distributed CB 
operation with Umited inter-base station communication, instantaneous information such as the channel 
vectors should not be exchanged since inter-base station communication delay is typically larger than 
the channel coherence time for mobile users. One possible way to roughly control the rate-tuple in the 
network is to design a table composed of sets of interference relaxation parameters, as in the right side 
of Fig. |7J based on the channel statistics. When the transmitters form a coordinating cluster, they can 
negotiate their rates based on the requests from their receivers for a communication session. In this phase, 
one set of interference relaxation levels from the table is picked, shared among the base stations, and 
used during the communication session. 

Heuristic guidelines to design the parameter table are based on the RZF parameterization itself. Note 
that ej = X^vVi'^^i ^'^ ® ^^ '■'^^ additional interference power relative to thermal noise power af at 
receiver i and e^ = 1 means that the SINR of receiver i is lower than the SNR of the same receiver 
by 3dB. Thus, the designed interference level should not be too high compared to the thermal noise 
level. Furthermore, to (roughly) obtain comer points of the Pareto boundary of the rate region, another 
heuristic idea works. One transmitter should use a nearly MF beam vector, and the rest of the transmitters 
should use nearly ZF beam vectors. More systematic ways based on vast computer simulation can be 
considered to design the parameter table. One possible way is as follows. We first generate a set of 
channel vectors randomly according to the channel statistics. For this realized channel set, we obtain 
graphs of interference relaxation parameters on the Pareto boundary. The process is repeated over many 
different channel reahzations and the best fitting graphs are obtained from the graphs of interference 
relaxation parameters of different channel realizations by some regression model. Finally, the table is 
constructed by selecting some points in the best fitting graphs. The parameter table in the right side of 
Fig- 13 is obtained in this manner for K = A^ = 2 when each element of channel vector is i.i.d. zero-mean 
complex Gaussian distributed with unit variance and the SNR is dB. Figures |7] (a) and (b) show the rate 
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Fig. 7. Performance of RZFCB with thie proposed rate control algorithms: The centralized rate control, marked with *, searches 
for the weighted sum rate maximizing point. (The weight vector w is shown in the figure.) The distributed rate control scheme, 
marked with +, sets the interference leakage levels as shown in the table. In Figs. |7] (a) and (b), 'virtual SINR' denotes the 
rate-tuple obtained by the virtual SINR (or SLNR) beamforming method in 1121 . 



control performance of the parameter table designed in this manner for two different channel realizations. 
It is seen that the heuristic method performs well; the five rate points are all near the Pareto boundary 
for each figure. 

Several advantages in the RZFCB are summarized below. 

• Real-time fully distributed operation is possible based on the proposed heuristic control approach. 
Transmitter i only needs to know {hji : j = 1, 2, • • • , K} and {cj^ : j = 1, 2, • • • , K}. 

• Once transmitter i knows {au, ■ ■ ■ , aj-i,j,ai-)-i,i, axi}. there exists a very fast algorithm, the SOPC 
algorithm, to design the RZFCB beam vector. Furthermore, in the case of K = 3, there is an approximate 
closed-form solution. 

• Transmitter i knows its SINR and achievable rate exactly, and its achievable rate is given by Ri = 

1 + /i+Vw^ • ^^' transmission based on this rate will be successful with high probability. This is 
true even when {ajj} are designed suboptimally, i.e., away from the Pareto boundary of the rate region. 
Thus, the RZFCB scheme is robust. 

• On the contrary to the ZF scheme, RZFCB does not require N > K. 



V. RZFCB FOR MIMO INTERFERENCE CHANNELS 

In this section, we consider the case that both transmitters and receivers are equipped with multiple 
antennas i.e., MIMO interference channels. In the MIMO case, we consider the weighted sum rate 
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maximization under the RZF framework and then propose a solution to the MIMO RZFCB based on the 
projected gradient method ||T61 . The rate control idea in the MISO case can be applied to the MIMO 
case too. 

A. Problem Formulation 

We assume that each receiver has M receive antennas and each transmitter has N transmit antennas. 
In this case, the received signal at receiver i is given by 

Yi = HiiVjSi + ^ YiifVjSj + rij, (30) 

where Hjj is the M x N channel matrix from transmitter j to receiver i, Vj is the N x di beamforming 
matrix, Sj is the di x 1 transmit symbol vector at transmitter i from a Gaussian codebook with Sj ~ 
CA/'(0, IrfJ, and iij ~ CJ\f{0,afl) is the additive noise. As in the MISO case, we have a transmit power 
constraint, ||Vj|||, < Pj, for transmitter j. The proposed RZF constraint in the MIMO case is given by 
an inequality with the Frobenius norm as 

||Hji Villi < Ojia], Vi, j / i (31) 

for some constant Oji > 0. As in the MISO case, the RZF constraints reduce to ZF constraints when 
Oji = for all i,j ^ i. With the MIMO RZF constraints, a cooperative beam design problem that 
maximizes the weighted sum rate is formulated as follows: 

Problem 4 (RZF cooperative beamforming problem): 

K 

1m + {afl + Bi)-iHiiV,Vf Hf , 

subject to (C.l) \\H.jiVi\\\ <ajia'j, Mi^j^i^ 

(C.2) \\V,\\l<Pu Vi, (32) 

where Wi > 0, J2i ^i = !> ^'^'^ ^« = SjV« ^ij^i^f^fj i^ ^^^ interference covariance matrix at 
receiver i. 

Note that, in Problem |4l the interference from other transmitters is incorporated in the rate formula 
through the interference covariance matrix Bj capturing the residual inter-cell interference under the 
RZF constraints. As in the MISO case, we will derive a lower bound on the rate of each user by 
exploiting the RZF constraints to convert the joint design problem into a set of separate design problems. 



max > Wi log 

{V.} ^^ 
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Note that, under the RZF constraints, the total power of interference from undesired transmitters is upper 
bounded as 



rl 



which impUes Bj < eiafl. 

Hassibi and Hochwald derived a lower bound on the ergodic rate of a MIMO channel with interference 
ll26l . However, their result is not directly applicable here since the rate here is for an instantaneous channel 
realization. Thus, we present a new lower bound under the RZF interference constraints in the following 
Lemma. 

Lemma 3: A lower bound on the rate of receiver i under the RZF constraints is given by 



log 



a 



1 



> log 



Im + 



T {l + e^ 



'-Hi ^ t ^ I '-'-II 



(34) 



where tr(Bi) = Y^j^i \\^ij^ jWf ^ ^i'^'i ^oi" ^^^ ^• 
Proof: The rate at receiver i is given by 



M 



log|I + *-^A,| = log J](l + Afc(*riAi)) 



(35) 



fc=i 



where *j = afl + Bj, Aj = HjjVjVf^H^, and Afc(X) denotes the fe-th largest eigenvalue of X. By 
the Rayleigh-Ritz theorem |[27l p. 176], we have 



^Mi^i^i) < -H^. 



X^$,X 



P^P 



<Ai(*-^A,) 



(36) 



for any non-zero vector x G C^^ and p := <^l x. From the Courant-Fischer theorem ll27l p. 179], the 
fe-th largest generalized eigenvalue of $t~^ Aj, A; = 1, • • • , M is given by 

p^*7^Aj*7^p 



Afc(*r'Aj 



max 
P^o, pec*', 
p±Pi,-- ,Pfc_i 



P^P 



(37) 



where pj is the eigenvector associated with the z-th largest eigenvalue of $^ ^-A-j$j ^ . Let Aj = 
UjSjUj^ be the eigen-decomposition of Aj, where 5]j = diag(Ai(Aj), • • • ,Am(Aj)). Then, for all 
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„ _ 1 _«. 

Afe(*riA.) = max ^ ^' 7^' ^ , 



p#o, pee^ P^P 

P^PI,-- ,pfc-l 



max „ 

P5^o, pec^^ P^P 

p±pi,-- ,pfc-i 



max 



z#0, zee", Z^Uf *jU,;Z 



J:.«/2tT I 

*■ UiZ_Lpi,--- ,pfc-i 



(z := Uf *, = p), 

(«) z-f^SiZ 

> max 



#0, ||z||=i, zec*^ z^Uf *,U,z' 



.if/2,, I 



E-=iA,(A,)|z,| 



max 



Z5^0, ||z|| = l, zGC", Z^Uf *,U,Z 

^fe + 1— 2fc + 2 — ■"— ^A/— 0, 

*■ UiZ±pi,--- ,Pfc-i, 
Ai(*,;) 

where (a) is satisfied since the feasible set for z is reduced and (b) is satisfied since ||z||^ = \z\\'^ + 

h \zk\^ = 1, Ai(Aj) > • • • > Afe(Aj), and z-^Uf ^jUiZ < Ai(*i) by Rayleigh-Ritz theorem. Based 

on dMl ). a lower bound on the rate is given by 

log|I + *-A.|>(l + MMV (39) 



Since *, = afl + E,y, Hi,V,Vf H,f , we have Ai($i) = af + Ai (^^^._,. H^.V^Vf H^^j , where the 
maximum eigenvalue of the interference covariance matrix is upper bounded by Ai ( Yljj^i ^ij Vj V^H^ j < 
trf Ylj^i ^ij^j^f^fj] ~ ^j^i ^ij'^f ~ ^i'^i- Thus, a lower bound of rate at receiver i is given by 

|I + *riAd > 



Im + 



<(l + e.)' 



Note that in (I34t the inter-user dependency is removed and the beam design can be performed at each 
transmitter in a distributed manner. Based on the lower bound ( [34l) . the RZFCB problem is now formulated 
as a distributed problem: 
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Fig. 8. Sum rate of RZFCB: (a) (K,M,N)=(3, 2, 6), (b) (K,M,N)=(3, 2, 8), and (c) (K,M,N)=(4, 2, 6). 



Problem 5 (The MIMO RZFCB problem): 

max (/>i(Vi) := log 
subject to (C.l) 

(C.2) 



Im + 2n I \ 






HjjVjII^ < ajiUj, 



Vj / i, 



^iWl < Pi, 



(40) 



for each transmitter i = 1, 2, • • • , K. 

Note that Problem|5]is now fully distributed. One of several known algorithms for constrained optimization 
can be used to solve Problem [5] for given {aji}. In particular, we choose to use the projected gradient 
method (PGM) by Goldstein f\M . The proposed PGM-based beam design algorithm for MIMO ICs is 
provided in Table Hill Detailed explanation of the beam design with PGM algorithm is provided in ||TT1. 

TABLE III 
Beam design algorithm for MIMO IC using PGM. 

For each transmitter i G {1, • • • , K}, 

0. Initialize Vj as the ZF beamforming matrix. 

1. Compute gradient of 0(Vj). 

2. Perform a steepest descent shift of V^. 

3. Perform successive metric projections of Vj onto constraint sets. 

4. Go to Step 1 and repeat until the relative difference of 0j(Vj) is 
less than a pre-determined threshold. 
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B. Numerical Results 

In this section, we provide some numerical results for the performance of RZFCB in the MIMO case. 
We consider three MIMO interference channels with system parameters {K, M, N) = (3, 2, 6), (3, 2, 8), 
and (4, 2, 6). In each case, we set aji = 0.01, 0.1 and 0.2 for all i and j. The step size parameter for the 
PGM is chosen to be 0.01 for all iterations. Figures |8](a), (b), and (c) show the sum rate performance 
of the ZFCB and RZFCB averaged over 30 independent channel realizations. In Fig. [8] (a) it is seen that 
the RZFCB outperforms the ZFCB at all SNR and the gain of the RZFCB over the ZFCB at low SNR 
is large when A^ = KM. This large gain at low SNR is especially important because most cell-edge 
receivers operate in the low SNR regime. In Fig. |8] (b) it is seen that the ZF scheme performs well when 
the number of TX antenna is more than enough and the dimension of ZF beams is large, as expected. 
In the case of A^ < KM as in Fig. (c), the ZFCB is infeasible but the RZFCB still works well. 

VI. Conclusion 

We have considered coordinated beamforming for MISO and MIMO interference channels under the 
RZF framework. In the MISO case, we have shown that the SOPC strategy with a set of well chosen 
interference relaxation levels is Pareto-optimal. We have provided (approximate) closed-form solutions 
for the SOPC strategy in the cases of two and three users and the SOPC algorithm in the general case for 
a given set of interference relaxation levels. In the MIMO case, we have formulated the RZFCB problem 
as a distributed optimization problem based on a newly derived rate lower bound and have provided an 
algorithm based on the PGM to solve the MIMO RZFCB beam design problem. We have also considered 
the rate control problem under the RZFCB framework and have provided a centralized approach and a 
fully-distributed heuristic approach to control the rate- tuple location roughly along the Pareto boundary 
of the achievable rate region. Numerical results validate the RZFCB paradigm. 
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